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Abstract 

Using the B s meson wave function extracted from non-leptonic B s decays, we reevaluate the 
rare decays B s — > £ + £~ 7, {£ = e, /i) in the standard model, including two kinds of contributions 
from magnetic-penguin operator with virtual and real photon. We find that the contributions from 
43 " magnetic-penguin operator 6 — ?■ S7 with real photon to the exclusive decays, which is regarded as 

qj to be negligible in previous literatures, are large, and the branchings of B s — > £ + £ 7 are nearly 

enhanced by a factor 2. With the predicted branching ratios at order of 10 , it is expected that 
the radiative dileptonic decays will be detected in the LHC-b and B factories in near future. 

PACS numbers: 12.15.Ji, 13.26.He 
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I. INTRODUCTION 



The standard model (SM) of electroweak interaction has been remarkably successful in 
describing physics below the Fermi scale and is in good agreement with the most experiment 
data. Thanks to the efforts of the B factories and LHC, the exploration of quark-flavor 
mixing is now entering a new interesting era. Measurements of rare B mesons decays such 
as B — > Xg'-f, B — > X s £ + £~(£ = e, //) and B s — > B s — > £ + £~^ are likely to provide 

sensitive test of the SM. In fact, these decays, induced by the flavor changing neutral currents 
which occur in the SM only at loop level, play an important role in testing higher order 
effects of SM and in searching for the physics beyond the SM 

n a. 

Nevertheless, these 

processes are also important in determining the parameters of the SM and some hadronic 
parameters in QCD, such as the CKM matrix elements, the meson decay constant /# s , 
providing information on heavy meson wave functions [sj]. 

The rare B inclusive radiative decays B — > Xg'-f and B — > X s £ + £~(£ = e,/i) as well 
as the exclusive decays B s — > have been studied extremely at the leading logarithm 

order 4] and high order in the SM and various new physics models. In previous works, 
prediction for the exclusive decays B s — > £ + £~~f have been carried out by using the light 
cone sum rule , the simple constituent quark model (6j|, and the B meson distribution 
amplitude extracted from non-leptonic B decays [7|. long distance QCD effects describing 
;he neutral vector- meson resonances <fi and J/^> family have received special attention in 
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10(. At parton level, B s — > £ + £ 7 decays have been thought to be obtained from decay 



b — > s£ + £^ , ~f, and further, from b — > s£ + £~ directly. To achieve this, a necessary work is 
attaching real photon to any charged internal and external lines in the Feynman diagrams 
of b —>■ s£ + £~ with two statements: i) Contributions from the attachment of photon to any 
charged internal propagator are regraded as to be strongly suppressed and can be neglected 
safely ; ii) Contributions from the attachment of real photon with magnetic- 

penguin vertex to any charged external lines are always neglected 

fla 

or stated to be 

negligibly small {3] • Here we would like to address that the conclusion of the first statement 
is correct, but the explanation is not as what it is described The second statement 

seems to be questionable, for that the pole of propagator of the charged line attached by 
photon may enhance the decay rate greatly which make some diagrams can not be neglected 
in the calculation. Since the weak radiative B-meson decay is well known to be a sensitive 
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probe of new physics, it is essential to calculate the Standard Model value of its branching 
ratio as precisely as possible. Although the second contribution has been calculated in Ref. 
8| , it mainly concentrated on the long distance effects of the meson resonances, whereas the 
short distance contribution which was incompletely analyzed. 

In this letter, we will concentrate on the short distance contribution to B s — > £ + £~7 
and check whether the contribution from magnetic-penguin operator with real photon to 
B s — > £ + £~ is negligible or not, and give some remarks including a comparison with other 
works. The paper is organized as follows. In sec. Ull we present the detailed calculation of 
exclusive decays B s — > £ + £~7, including full contribution from magnetic-penguin operator 
with real photon. Sec. IIHI contains the numerical results and comparison with previous 
works, and the conclusions are given in sec. IIVI 

II. THE CALCULATION 

In order to simplify the decay amplitude for B s — > we have to utilize the B s meson 

wave function, which is not known from the first principal. Fortunately, many studies on 
13] and B s decays [l^] have constrained the wave function strictly. It 



non-leptonic B 
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was found that the wave function has form 

§b s = (0b s + m Bs )^ 5 <t> Ba (x), (1) 

where the distribution amplitude 4> Ba (x) can be expressed as flol ]: 

<p Ba (x) = N Ba x\l -x) 2 exp (2) 

with x being the momentum fractions shared by s quark in B s meson. The normalization 
constant N Ba can be determined by comparing 

(0 \sj^ 5 b\ B s ) = i / cf) Ba (x)dxTr [7 M 7 5 (^s s +m Ba )j 5 ]dx = -Aip B I (f) Ba (x)dx (3) 
Jo Jo 

with 

{0\sn 5 b\B s ) = -if Ba p Bs , (4) 
the B meson decay constant f Ba is thus determined by the condition 



f l 1 

^ 0U s (x)rfx = -f Bs . 



(5) 



Let us start with the quark level processes b — > s£ + £~. They are subject to the QCD 
corrected effective weak Hamiltonian, obtained by integrating out heavy particles, i.e., top 
quark, higgs, and W ± , Z bosons: 



U ef f(b-> s£ + £~) 



OtemGF 



7T 



v tb v* 



2Cf f m b _ 



sia» u q u P R b + C e 9 ff sYP L b 



(6) 



where Pl,r = (1 =F 7s)/2, q 2 is the dilepton invariant mass squared. The QCD corrected 
Wilson coefficients C% , Cg and C\q at /i = m;, scale can be found in Ref. 16]. 

If an additional photon line is attached to any of the charged lines in diagrams con- 
tributing to the Hamiltonian above, we will have the radiative leptonic decays b — > s£ + £~^. 
Therefore, there are two kinds of diagrams: photon connecting to the internal propagators, 
and photon connecting to the external line. As addressed in the introduction, the contribu- 
tion from the first kind of diagrams is neglected safely. Now we will only consider the second 
category of diagrams which are displayed in Fig. [TJ At first, we recalculate the diagrams 
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FIG. 1: Feynman diagrams for B q ->• £ + £~-f (q = d, s) in the SM . The black dot in (a), (c) and (e), 
(f) stands for the magnetic-penguin operator O7 with virtual and real photon, respectively, and 
the black dot in (b),(d),(g) and (f) denotes operators O9 and 0\q. 



(a)-(d) in Fig. [T]with photon emitted from the external quark lines b or s by using the B s 
meson wave function extracted from non-leptonic B s decays. Note that these diagrams are 
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already studied in previous literatures and considered as giving the dominant contribution 
to B s — > At parton level, the amplitudes for transition b — > si + i~^j can be calculated 

directly using from the Hamiltonian of b — > s£ + £~ in For example, contribution from 
the magnetic-penguin operator with virtual photon shown in Fig. [T] (a) reads: 

A a = -iGee^C^ (s [ 2 ^ + /V &) [^] , (7) 

where p& iS , fc denotes the momentum of quarks and photon respectively, e is the vector 
polarization of photon and G = a em GFVtbV t */(\/27i), e d = —1/3 is the number of electrical 
charge of the external quarks. In deriving above equation, we have used motion equation 
for quarks and q^i^^i = 0. Using Eqs. ([!]) and ([7]), we write the amplitude of B s — > £ + £~7 
at meson level as: 

A a = 2zGee d ^^CV e// ^- f ^-dx [k,q -e-e.k-q- ie^fcV] [If 1] ,(8) 
q PB 3 -kJ x 

where x, y = 1 — x are the momentum fractions shared by s, b quark in B s . By doing the 

similar calculation to diagrams (b)-(d), the decay amplitude is then obtained as: 



A a+b+c+d = iGee d —^- - j \C 1 ie a ^ v p° B e ti k lJ + C 2 p Ba (e IM k u - k^e u ) 

+C W [C + ie a ^ v p B /k v + C.p Bs (e^k u - k„e v )] (9) 
The form factors in Eq. (Q are found to be: 



C eff _ 2 H^ C eff 



C 2 = Cf'CL - 2 1 ^^C e 7 ff C + , (10) 

where 



q 2 



C± = J (^±~\<f> Bs (x)dx. (11) 

The expression in fl9]) can be compared with Ref. Q. 

Now we will focus attention on calculating the diagrams (e) and (f) in Fig. [1] which are 
always neglected in other works. In these two diagrams, photon of the magnetic-penguin 
operator is real, thus its contribution to B s — > is different from that of magnetic- 

penguin operator with virtual photon in diagram (a) and (c). We get the amplitude: 

A e+f = i2Gee d C e 7 ff mbm 2 Bs 1 C + [k^q -e-e^k-q- ie, va ^k a q?] [l-fl] , (12) 

Q Pb 3 'Q 



with coefficients C + obtained by a replacement: 
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(13) 



where z 



2pB s q 



. Note that pole in C + corresponds to the pole of the quark propagator 



when it is connected by the off-shell photon propagator. Thus the C + term may enhance 
the decay rate of B s — > and its analytic expression reads 



C, 



2N Bs Triz 2 (l - z) 2 exp 
1 
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(14) 



As the contribution from the FigJT] (g) and (h) with photon attached to external lepton 
lines, considering the fact that (i) being a pseudoscalar meson, B s meson can only decay 
through axial current, so the magnetic penguin operator O-j 's contribution vanishes; (ii) the 
contribution from operators Og, 0\q has the helicity suppression factor m^/m Ba , so for light 
lepton electron and muon, we can neglect their contribution safely. 

The total matrix element for the decay B s — > is obtained a sum of the A a+ b +c+ d 

and A e+ f. After summing over the spins of leptons and polarization of the photon, and 
then performing the phase space integration over one of the two Dalitz variables, we get the 
differential decay width versus the photon energy E 7 , 
dT 



01 F It/ t/ *,2 



108tH 



V tb V t *f(m Bs - 2E 1 )E 1 [\C,\ 2 + |C 2 | 2 + C 2 (|C + | 2 + |C_| 2 )] . (15) 



The coefficients C* (i = 1, 2) can be obtained by a shift: 

2m b m Bs p Bs ■ k^- 



Pb s ■ q 



C+C e 7 ff . 



(16) 



III. RESULTS AND DISCUSSIONS 



The decay branching ratios can be easily obtained by integrating over photon energy. In 
numerical calculations, we use the following parameters 181 ] : 

1 



a,. 



137 



, G F = 1.166 x 10~ 5 GeV" 2 , 



m b 



4.2GeV, 



\V tb \ = 0.88, \V ts \= 0.0387, \V td \ = 0.0084 

m Bs = 5.37GeV, co Ba = 0.5, f Ba = 0.24GeV, r Bs = 1.47 x 10~ 12 s. 

m B o = 5.28GeV, cj Bd = 0.4, f Bd = 0.19GeV, r Bd = 1.53 x 10~ 12 s. 

The ratios of B s — > £ + £~j with and without the contribution from magnetic-penguin oper- 
ator with real photon are shown in Table [I] together with results of B^ s — > i + £~ r y from this 
work and other models for comparison. The errors shown in the Table [I] comes from the 
leavy meson wave function, by varying the parameter u Bd = 0.4 ± 0.1, and u Bs = 0.5 ± 0.1 
3]. Note that, the predicted branching ratios receive errors from many parameters, such as 
meson decay constant, meson and quark masses etc. 



TABLE I: Comparison of branching ratios with other model calculations 



Branching Ratios (xlO 9 ) 


Our Results 


Quark Model 


light cone 




Excluded Fig.l(e),(f) 


Included Fig.l(e),(f) 


Ref.[7] 


Rcf.[6] 


Ref.yy 


B s e+£~j 


o 7 4 +1.76 
— 1.00 

n i6 +an 

u - 1D -0.05 


7.45^;1 
0.31^ 


1.90 
0.08 


6.20 
0.82 


2.35 
0.15 



A couple of remarks on the B s rare exclusive radiative decays are follows: 

1. As pointed out in Ref. , the branching ratios are proportional to the heavy 
meson wave function squared, the radiative leptonic decays are very sensitive probes 
in extracting the heavy meson wave functions; 

2. The contributions from magnetic-penguin operator with real photon to the exclusive 
decay is large, and the branching of B s — > £ + £~ , ~f is enhanced nearly by a factor 2 
compared with that only contribution from magnetic-penguin operator with virtual 
photon and nearly up to 10 -8 , implying the search of B s — > £ + £~ r j can be achieved in 
near future. 

3. Due to the large contributions from magnetic-penguin operator with real photon, the 
form factors for matrix elements (7|s7 /x (l ± 7s)6|i? s ) and i^j\sa^ v {l ± r y? ) )q u b\B s ) as 
a function of dilepton mass squared q 2 are not as simple as l/{q 2 — q^) 2 where q 2 
is constant jl7]. The B s — > 7 transition form factors predicted in this works have 
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also some differences from those in Ref. js-10]. For instance, Ref. predicted the 
form factors -Fry(g 2 ,0), Fta(q 2 ,0) induced by tensor and pseudotensor currents with 
emission of the virtual photon, as shown in diagrams (a) and (c) of FIG. [U are only 
equal at maximum photon energy, whereas the corresponding formula in this work 
have the same expression as — ^-^-C + oc l/(q 2 — qfi) in Eq. OH]). The research of 
B s — > £ + £~'-f may give some hints on these form factors. 

At this stage, we think it is necessary to present a few more comments about the calculation 



of Ref. js]. In order to estimate the contribution of emission of the real photon from the 
magnetic-penguin operator, the authors of Ref. jg] calculated the form factors Fta,tv(0, q 2 ) 
by including the short distance contribution in q 2 — > limit and additional long distance 
contribution from the resonances of vector mesons such as p°, u for Bj decay and (f> for 
B s decay. Obviously, this means the pole mass enhancement of the valence quark were not 
appropriately taken into account. Moreover, if -Ftm,tv(0, q 2 ) = -Pta.tv (0,0) stands for the 
short distance contribution, it seems double counting since in this case photons which emits 
from magnetic-penguin vertex and quark lines directly are not able to be distinguished. 

IV. CONCLUSION 

We evaluated the rare decays B s — > / ~y£ + £~ in the SM, including two kinds of contributions 
from magnetic-penguin operator with virtual and real photon. In contrast to the previous 
works which treated contribution from magnetic-penguins operators with real photon to the 
decays as negligible small, we found that the contributions is large, leading to the branching 
of B s — > £ + £~ r y being nearly enhanced by a factor 2. In the current early phase of the 
LHC era, the exclusive modes with muons in the final states are among the most promising 



decays. The decay B s — > is likely to be confirmed before the end of 2012 [19|. Although 

there are some theoretical challenges in calculation of the hadronic form factors and non- 
factorable corrections, with the predicted branching ratios at order of 10~ 8 , B s — > £ + £~~ / -f can 
be expected as the next goal once B s — > n + fj>~ measurement is finished since the final states 
can be identified easily and branching ratios are large. Our predictions for such processes 
can be tested in the LHC-b and B factories in near future. 
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